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Dependence of Modulation Amplitude on Electron Density in Unidirectional Lateral 
Superlattices: The Effect of the Thickness of the Two-dimensional Electron Gas 

Akira Endo * and Yasuhiro I ye 
Institute for Solid State Physics, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8581 

The amplitude Vo of unidirectional periodic potential modulation introduced by a surface 
grating into a two-dimensional electron gas (2DEG) formed at AlGaAs/GaAs heterointerface 
is measured as a function of electron density n e by analyzing commensurability oscillation of 
the magnetoresistance. The electron density is varied either by applying a bias to a metallic 
back gate or by illumination. The amplitude decreases with increasing density, with the rate 
\dVo /dn e | roughly an order of magnitude larger for the former method. The result is interpreted 
in terms of the rate, dEi/d(5E c ), of the change in the first subband level E\ in response to 
the variation of the conduction-band edge SE C above the heterointerface. The rate crucially 
depends on the thickness of the 2DEG. 

KEYWORDS: two-dimensional electron gas, AlGaAs/GaAs single heterostructure, potential modulation, 
lateral superlattice, electron density, thickness, back gate, illumination, commensurability 
oscillation 



1. Introduction 

A two-dimensional electron gas (2DEG) formed at 
Al^Gai-^As/GaAs heterointerface has been the basis of 
a multitude of experimental studies for low-dimensional 
electron systems. The devices for such experiments are 
produced by processing the surface of the 2DEG wafers 
by micro- (or nano-) fabrication techniques: depositing 
patterned gates or etching off part of the wafers, thereby 
introducing electrostatic potential modulation or con- 
fining electrons to smaller areas. 1 The modern nano- 
fabrication technology grants fine control of the lateral 
dimensions of the patterns down to the length scale less 
than 100 nm, limited only by the depth of 2DEG plane 
from the surface. In contrast to the high controllabil- 
ity of the lateral length scale, the energy scale involved, 
the magnitude of the potential modulation, allows much 
poorer management. The amplitudes are often adjusted, 
chiefly on empirical basis, only for a certain particular 
purpose, e.g., depleting electrons underneath the metal- 
lic gate by applying negative bias. It is usually not easy 
even to measure, much less to precisely adjust, the am- 
plitude of the modulation for generic cases. An impor- 
tant exception is when the modulation possesses a uni- 
directional periodicity. The unidirectional lateral super- 
lattice (ULSL) reveals the amplitude Vq of the modula- 
tion through the amplitude of the commensurability os- 
cillation (CO), the magnetoresistance oscillation result- 
ing from the commensurability between the period a of 
the modulation and the cyclotron radius R c . 2 

In the present paper, we report the behavior of the 
modulation amplitude when the electron density n e is 
varied, in pursuit of better quantitative understanding of 
the energy scale relevant in the study of low-dimensional 
electron systems. The amplitude Vo is measured using 
CO of ULSL. We shed light on the effect of finite thick- 
ness of the 2DEG wave function, which is often neglected 
in the first order approximation but nevertheless can be 



of vital importance for phenomena in which electron- 
electron interaction plays a crucial role. 3 In a ULSL, 
a grating that introduces potential modulation to the 
2DEG plane is placed on the front surface. Therefore, it 
is clear that a metallic uniform front gate, a standard 
device to vary n e , cannot be employed without heav- 
ily affecting the modulation amplitude. We have instead 
made use of metallic gate on the backside of the wafer, 
or illumination by light emitting diode (LED). Naively, 
one would expect the modulation amplitude to be basi- 
cally insensitive to the back gate, since the grating and 
the back gate, being situated on the opposite side of the 
2DEG, is expected to cause minimal crosstalk. It turns 
out, however, that the amplitude varies substantially in 
response to the back gate, much more pronounced than 
when n e is modified by illumination. The large variation 
is ascribed to the modification in the width of the well 
confining the electrons in the vertical direction, hence to 
the thickness of the 2DEG wave function. The experi- 
mental methods and results are described in §2 and 3, 
respectively. The interpretation of the experimental re- 
sults is discussed in detail in §4, followed by concluding 
remarks in §5. 

2. Experimental 

Two ULSL samples with differing periods were pre- 
pared from the same Al^Gai-^As/GaAs (#=0.33) single- 
heterostructure (SH) 2DEG wafer having the mobility 
/x~70 m 2 /Vs and electron density n e — 2xl0 15 m~ 2 at 
T—4.2 K. The structure of the wafer was (from the 
front surface) 10 nm GaAs cap layer, 40 nm Si-doped 
(A^ S i=2xl0 24 m~ 3 ) ALGai_ x As layer, 40 nm undoped 
Al x Gai_ x As spacer layer, and 1 /im GaAs layer with 
2DEG channel residing near the interface to the upper 
layer. As depicted in Fig. 1(c), a pair of Hall bars, 44x 16 
/im 2 , were defined in series by wet-etching; one of these 
was further processed into ULSL and the other was re- 
served as reference. The potential modulation was intro- 
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Fig. 1. (a) A typical magnctoresistance trace showing PMR 
and CO (taken at 4.2 K). (The example is for sample B at 
n e =2.2xl0 15 m~ 2 ). Arrows mark the position of -Bj n f. Minima 
in CO are labeled by their index numbers, (b) Oscillatory part of 
CO, obtained by subtracting the slowly varying background. The 
dotted curve represents a fit to cq. (1). (c) Schematic drawing of 
the sample. 



duced employing strain-induced piezoelectric effect 4 by a 
grating of high-resolution negative electron-beam resist 5 
placed on the front surface. 6 The periods are a=161 and 
138 nm for samples A and B, respectively. A uniform 
metallic (gold) gate was deposited on the rear side of 
the wafer after thinning the wafer from 0.3 mm down 
to about 0.1 mm by wet-etching in order to enhance 
the effectiveness of the back gate. No signs of deterio- 
ration in the mobility were observed after these process- 
ing. 7 The electron density varies, by Anb g , in response 
to the bias V bg applied to the back gate. From the rate 
d(An bg )/dVb g =5.0xlO- 3 and 4.1xl0~ 3 in 10 15 m^/V, 
the distance of 2DEG plane from the back gate was es- 
timated to be (ibg=140 and 170 /xm by a simple capaci- 
tance model eAnb g =(eoe/(ibg)Vb g for samples A and B, 
respectively. 

Figure 1(a) shows a typical magnetoresistance trace 
of our ULSL, taken by standard low-frequency ac lock-in 
technique at 4.2 K. The trace exhibits CO as well as pos- 
itive magnetoresistance 8 (PMR) emanating from B=0. 
The amplitude Vb of the unidirectional periodic poten- 
tial modulation 9 V mo d{x) = Vb cos(qx) with q = 2ir/a, 
can be deduced from the amplitude of CO. In order to 
obtain accurate quantitative value, however, care should 
be taken of the impurity scattering that scatters elec- 
trons out of the cyclotron orbit before completing a cycle. 
Such scattering reduces the CO amplitude, more effec- 



tively for smaller B where cyclotron circumference 2irR c 
is large, and can lead to the underestimation of Vq if the 
standard prescription 10 ' 11 is used for the analysis. The 
present authors have shown 6 that the scattering can be 
well accounted for by the inclusion of an additional factor 
A(n/fiwB) in the formula as, 

Ap° x S x C / 7T \ / 

P0 ~ Uw-B/ \Ta 



1 



1 m 2 Vi 



IB I sin 2tt 



2R C 



where A(x)=x/ sinh(x), kBT a ~(l/2n 2 )(akp /2)Huj c with 
oj c =eB/m*, $ =Ve, ^=eh / 2m* , and R c =hkp/e\B\ 
with fci?=(27rn e ) 1 / 2 the Fermi wave number. As shown in 
Fig. 1(b), the oscillatory part of the magnetoresistance 
can be fitted, with and Vb as fitting parameters, to 
eq. (1) very well. The parameter nw describes the degree 
of decay of CO amplitude with decreasing B, and was 
shown to be able to be identified with 6 the single-particle 
or quantum mobility /iq deduced from the analysis of the 
decay of Shubnikov-de Haas oscillation 12 in the adjacent 
reference plain Hall bar. An alternative way to deduced 
Vb is provided by the analysis of PMR. The present au- 
thors have recently pointed out 13 that inflection points 
B[ n { at which curvature of the magnetoresistance changes 
from concave down to concave up, as pointed by the ar- 
rows in Fig. 1(a), can be found on the PMR when Vb 
is small enough, and B[ n f corresponds to the extinction 
field, 



Br. 



2im*V 



(2) 



where streaming orbits, the orbit of electron confined in a 
single period, cease to exist. The values of Vb deduced by 
this second method agree very well with those drawn us- 
ing eq. (1), confirming that reliable values are obtained. 

3. Results 

Figure 2(a) shows Vq obtained using eq. (1) from the 
amplitude of CO for samples A and B, plotted as a 
function of n e . The electron density n e is varied either 
by back-gate voltages (solid symbols) or by illumination 
with an infrared LED resorting to the persistent photo- 
conductivity effect (open symbols). The ratio of Vb to 
the Fermi energy Ep ranges from 2 to 4% for sample A, 
and 1 to 2% for sample B. It can clearly be seen that Vb 
decreases with increasing n e for both samples. The rate 
of the decrease dVo / dn e is nearly an order of magnitude 
larger when n e is varied by the back gate. Although Vb of 
sample A is much larger than that of sample B, 14 the rel- 
ative change, Vo normalized by the value at n e =2.0x 10 15 
m~ 2 , shows almost the same behavior for both samples, 
as demonstrated in Fig. 2(b). 

Decrease of the modulation amplitude V m response 
to the increase in n e by a back gate was also reported by 
Soibel et at 15 As will be detailed in the next section, the 
decrease cannot be ascribed to the screening by 2DEG, 
which is independent of n e unless a is less than the half 
of the Fermi wavelength \p. In the case of rcf. 15, Vb, 
hence the ratio Vo/Ep, is roughly an order of magni- 
tude larger than those of our samples. The authors of 
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Fig. 2. (a) Amplitudes Vo of potential modulation, deduced from 
CO at 4.2K using eq. (1), plotted as a function of the electron 
density n e . Solid and open symbols indicate that the densities 
were varied by the back gate or by illumination, respectively. Up- 
ward and downward triangles represent samples A (a=161 nm) 
and B (a=138 nm), respectively, (b) Amplitudes normalized by 
interpolated (n e back-gate controlled) or extrapolated (n e illu- 
mination controlled) values at n e =2.0 X 10 15 m~ 2 . Calculated 
curves simulating the variation of n e by the back gate and il- 
lumination arc shown by solid and dotted curves, respectively, 
for A r J 4=l-0 and 2.0xl0 20 m — 3 as thin and thick curves, respec- 
tively. (See text for details.) 



ref. 15 pointed out that their Vo is not small enough to 
be treated within the standard framework of the screen- 
ing theory; they claimed that alteration of the 2DEG 
density of state by their large Vo should be taken into 
consideration. Obviously, such explanation does not ap- 
ply to our samples having Vo of only a few percent of Ep. 
In the next section, we will describe our interpretation 
of the behavior of Vo caused by the change of n e . 

4. Discussion 

4-1 Screening 

We start by examining the effect of screening. In 
the Lindhard formula, the dielectric constant for an 
ideal zero-thickness 2DEG is given by 16,17 e ee s (<z) with 
e =8.85x 10~ 12 F/m the permittivity of the vacuum, 
e the relative dielectric constant of the host crystal 
(e=13.18 for GaAs), and 



e s (q) = 1 + 



a B q 



1 - Re 



2k F 

q 



(3) 



expressing the screening by the 2DEG. 
a* B =4:TT€oefi 2 /m*e 2 represents the effective Bohr ra- 
dius for the host material (10.4 nm for GaAs). In eq. 
(3), only the last term depends on n e , which disappears 
at q<2k F . Namely, the change in n e , per se, does not 
affect the screening for a>X F /2, where e s (q) equals 
e TF(<?) = l + 2/a|j<7 calculated by Thomas-Fermi ap- 
proximation. For the range n e =1.5-3.0x 10 15 m ~ 2 of 
the electron density investigated in the present study, 
Xp =45-65 nm, and therefore, the screening by 2DEG 
of our potential modulation V mo d(x) will not depend 
on n e . For realistic 2DEG with finite thickness, £tf(?) 
should be slightly modified to include the form factor 
F(q): 



tTF(q) = l + - r -F(q). 

a B q 

The form factor is defined as 18-20 



(4) 



F(q) 



dk f(k)f(-k) 
2tt q 2 + k 2 



2 "I 

If dzdz'f(z)f(z')aw(-<l\z - z% (5) 



using the Fourier transform f(k) of the (normalized) 
electron distribution function f{z)=\^{z)\ 2 in the z 
direction, the direction normal to the 2DEG plane. 
ip(z) designates the envelope function of the electron 
wave function ip(z) exp{i(k x x + k y y)\. It can readily 
be seen that F(q) = l for f(z)=S(z), an ideal 2DEG. 
As a rule of thumb, F(q) decreases when f(z) is 
widely spread: the screening is less effective for thicker 
2DEG. For example, for Fang-Howard wave function 16 ' 21 
V'FH(2)=e(z)(6 3 /2) 1 / 2 zcxp(-6z/2), F(q)=b(8b 2 + 9bq + 
3q 2 )/8(b + q) 3 = l-W/8(q/b) + 3(q/b) 2 + 0[(q/b) 3 ].RcTC, 
O(z) represents the unit step function. The variational 
parameter b, given by &=[487r(nd cp i + lln e /32)/a| ? ] 1 / 3 
with ndepi the depletion charge, is inversely proportional 
to the thickness of the wave function [rms thickness 
((z 2 ) - (z) 2 ) 1 ' /2 =3 1 ' /2 /6], and therefore F(q) decreases 
with increasing thickness. 

The thickness of 2DEG varies with n ei and therefore 
the change in n e can, in principle, alter the modulation 
amplitude seen by the electrons through the change in 
the screening. However, it is obvious that this thickness- 
mediated change in the screening is unable to explain 
our Vo — n e relationship even in the qualitative level. As 
we will see in the next section (Fig. 5(a)), the rms thick- 
ness decreases with increasing n e , leading the screening 
to be more effective thus letting V) smaller, if n e is var- 
ied by illumination (or by a front gate). However, it is 
the other way around when the back gate is used: the 
thickness, hence Vo, is expected to increase with increas- 
ing n e , at variance with the experiment. Furthermore, 
it can be demonstrated that the change in £tf(<z) is too 
small, roughly 1% at most, by evaluating eq. (5) using 
numerically calculated ">p(z) that will be described in the 
next subsection. The observed Vb vs n e requires a mech- 
anism other than the screening for its explanation. In 
what follows, we neglect altogether the change in €TF{q) 
by n e . 
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4-2 The conduction-band edge and the subband level 

For the range of n e and the temperature, 4.2 K, con- 
sidered in the present paper, electrons occupy only the 
lowest subband in the confinement potential. The po- 
tential modulation seen by the electrons is the spatial 
variation of this subband level, E\, with respect to the 
Fermi level. By contrast, it is the conduction-band edge 
E c , not Ei, that is modified by the external devices to 
introduce potential modulation (a grating placed on the 
surface in the present study); they do not directly cou- 
ple to Ei. The distinction is not important when the 
two energy levels shift in parallel, which will be approxi- 
mately the case for 2DEG in a narrow quantum well. In 
general, however, the shift in Ei does not necessarily ex- 
actly follow the shift in E c . The subband level Ei , hence 
the conversion rate dEi/d(5E c ) from the shift 8E C in the 
conduction-band edge, evaluated just above the hetcroin- 
terface where wave function of the electron vanishes, to 
Ei , is dependent on the width or the profile of the con- 
finement potential V(z) in which the subband is formed. 
In ALrGai-^As/GaAs SH 2DEG, electrons themselves 
and the back-gate voltages are important ingredients for 
determining the profile of V(z), as we will see below. 
Therefore, dEi/d(SE c ) varies with electron density or 
back-gate conditions. In what follows, we demonstrate 
by numerical simulation using simplified model that the 
observed variation of Vb by the back gate or illumina- 
tion is mainly attributable to the concomitant change in 
dEi/d(5E c ). 

The potential well confining the electrons is given 

by 22,23 



and Posisson's equation, 



V(z) = V b e(-z) + (-e)<f>(z) + V xc (z). 



(6) 



Here, we define z=Q as the heterointerface, and the z-axis 
points down into the substrate (z<0: Al^Gai-^As and 
z>0: GaAs). V(z) is comprised of the band discontinuity 
V b Q(-z) at the heterointerface (14=292 meV for £=0.33 
at 4.2 K) and the electrostatic potential (— e)<f>(z) origi- 
nating from the ionized donor (Si) at doped Al^Gai-^As 
layer, remnant background impurities, as well as the 
electrons themselves (the Hartree term). The exchange- 
correlation effects, the effects of electron-electron inter- 
action beyond the Hartree approximation, is also taken 
into account as parameterized potential 24 V xc (z) after 
Stern and Sarma, 23 



V xc (z) 



r s ( 21 \ 
1 + 0.7734- In (l + -J 



(4/9 7 r) 1 /3 7rrj 



where 



r s = r s {z) 



^ira* B 3 n e \tp(z)\ 2 



-1/3 



-V, 

(7) 



(8) 



and Ry*=e 2 /87reoe«|j the effective Rydberg energy (5.2 
meV for GaAs) . The electrostatic potential is determined 
by solving the Schrodingcr equation, 



dz 2 ' 



(~e) 



[n e \i>{z)\ 2 + N A Q(z)} (10) 



self-consistently. Here we have made several types of sim- 
plifications: the difference in the effective mass and di- 
electric constant between the two materials is neglected 
and those for GaAs (TO*=0.067m e and e=13.18) are used 
throughout, and accordingly the image potential energy 
is ignored; possible atomic-scale smooth grading of the 
band discontinuity at the interface considered in ref. 23 is 
neglected; the background residual impurity Na is taken 
into account only at z>0. 25 These approximations are 
expected not to affect the present argument very much. 
In GaAs-based materials grown in the modern molecu- 
lar beam epitaxy (MBE) machines, it is known that the 
residual impurity is mainly composed of carbon that is 
unintentionally incorporated into the crystal. In GaAs, 
the carbon works as an acceptor, whose density Na in the 
particular 2DEG wafer used in the present study is diffi- 
cult to know. It can be estimated, however, to be close to 
1.7xl0 20 m~ 3 inferred from Hall measurement for nom- 
inally undoped GaAs bulk crystal grown slightly before 
the present 2DEG wafer was grown, using the identical 
MBE chamber. Therefore we use Na near this value for 
calculations. 

The effects of the ionized donor, the background im- 
purity, and the back gate are taken into account as 
suitable boundary conditions. The slope of the electro- 
static potential just below the 2DEG, i.e., where |f/'(- 2; )| 2 
vanishes at (z)<Z"CZdep] (with (z) and z^ cp i denoting 
the average position of the wave function and the de- 
pletion layer thickness, respectively), is determined by 
the electric field due to the depletion charge ridcpi and 
the slope by the back gate. We parameterize the ef- 
fect of the back gate by the change in the electron 
density Arthg induced by applying the back-gate volt- 
age Vb g ; eAnb g =(eoe/dbg)14g as defined before. Thus, 
by Gauss's theorem, — d(j>/dz={— e/eoe)(n<iepi — Arib g ). 
The slope just above the 2DEG (where |t/>(z)| 2 vanishes 
at — d sp <^z<0, with d sp =40 nm representing the thick- 
ness of the spacer layer) results from the charges from 
the donor (or equivalently from the charges from 2DEG 
electrons in addition to the depletion charge and back- 
gate contribution, considering the charge neutrality), and 
therefore — dcp/dz=(— e/eoe)(ndcpi + n e — Anb g ). The de- 
pletion charge and the depletion layer thickness are found 
by the conditions (f>(zdepi)=E g /(— e) and d(f>/dz(zd cp \)=0, 
setting Anb g =0. We obtain ndc P i=(2eoe./V aE 9 j e 2 ) 1 / 2 and 
Zdcpi—{^o^E g /e 2 NA) 1 ' 2 , respectively, the latter being 
the order of micro meters. As a band gap of the GaAs, 
we adopted E g =1.52 eV, the value at 4.2 K. The solu- 
tion of eqs. (6), (9), and (10) with the above boundary 
conditions are formally written as, 



<t>(z) 



(n dcp i - An hg )z 



-e(z)±N A z 2 + n e {I(z) - I(z dcp i)] j , (11) 




z(nm) z(nm) 



Fig. 3. Numerically calculated conduction-band edge E c and the 
wave function ip, for iV J 4=2xl0 20 m — 3 . The heterointerface re- 
sides at z=0. Calculations were done for the values of n e (in 10 15 
m — 2 ) from 1.5 to 3.0 by the increment of 0.1 each, and displayed 
with progressively brighter grayscale. The back gate was set (a) 
neutral (An bg (in 10 15 m~ 2 )=0.0) (b) positive (+0.2), and (c) 
negative (—0.5). 



with 

I(z) = z- f dz' f ' dz"\^P{z")\ 2 . (12) 

J — oo J — oo 

Numerical solutions of the wave function tjj(z), and the 
z dependence of the conduction-band edge E c , which is 
no other than the confinement potential V(z), is plotted 
in Fig. 3. Figure 3(a) to (c) correspond to different set- 
tings of the back-gate voltages Anb g =(eoe/edb g )Vb g . For 
each Anbg, calculations are repeated for different n e 's 
ranging from 1.5 to 3.0xl0 15 m~ 2 with an increment of 
O.lxlO 15 m~ 2 each. As a matter of course, the confine- 
ment potential, hence the wave function, becomes thicker 
(thinner) for a positive (negative) back-gate voltage. For 



Fig. 4. The conduction-band edge E c near the heterointerface 
and the lowest-subband level E\. The values of Na, n e , and 
An^, an d the grayscale are the same as in Fig. 3 

a fixed Aribg, the width of ip(z) decreases with increasing 
n e , owing to the last term in eq. (11). The change in the 
thickness is more conspicuous for larger (more positive) 
Anbg- This is more clearly illustrated in Fig. 5(a), which 
plots rms thickness of the wave function versus n e for 
varying settings of Arib g ranging from —0.5 to +0.2 with 
an interval of 0.1 (in 10 15 m~ 2 ). Figure 4 shows close-up 
of E c in the proximity of the heterointerface along with 
the lowest subband level E\. The conduction-band edge 
E c at the heterointerface, z=0, displays downward shift 
with the increase of n e . Accompanying this deepening of 
the bottom of the confinement potential, the subband 
level Ei is also displaced downward but with smaller 
decrement. Close inspection of Fig. 4(a)-(c) reveals that 
the downward shift in E c by n e is larger (smaller) for 
more positive (negative) back-gate voltages, while the 
extent of shift in E\ remains virtually unchanged. E c at 
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Fig. 5. (a) rms thickness of the wave function vs n e , (b) the shift 
AEi of the lowest subband level E\ (measured from the value 
at n e =2.0xl0 15 m~ 2 ) vs the shift 8E C of the conduction-band 
edge just above the heterointerface, (c) the rate of the change of 
Ei in response to the change in SE C , for different settings of the 
back gate (Anj,g from —0.5 to +0.2xl0 15 m~ 2 at the interval 
of 0.1 xlO 15 m -2 , each plotted by differently shaped symbols), 
and 7V J 4=2xl0 20 m -3 . The large (small) solid circles designate 
how the values evolve when n e is varied by the back gate (by 
illumination). The right axis in (c) shows the ratio with the value 
at n e =2.0x!0 15 m" 2 and An bg =0.0x 10 15 m" 2 . 



the heterointerface approximately represents the E c just 
above 2DEG, the E c that we will be looking at (to be de- 
noted as SE C ), apart from minor modification caused by 
not including the contribution from the tail of the wave 
function in the Al^Gai-^As barrier layer and including 
the contribution from V xc (0). More quantitative account 
is given in Fig. 5(b). 

Before explaining Fig. 5 further, we look into more de- 
tail of the behavior of the last term in eq. (11), which gov- 
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cms the n e dependences we have just pointed out. The 
function I(z) appearing in the term tends to z for z<C0, 
and to a constant value dependent on the distribution of 
ip(z) for large enough z(~^>(z}), as exemplified in Fig. 6 
for particular values of N^, n e , and Arib g - It can easily 
be perceived from the definition eq. (12) that the con- 
stant value at z~^>(z), 7(oo)[~/(zdepi) since Zdcpi^M 2 )], 
increases with increasing width of ip(z). Physically, this 
reflects the property of the solution of the Poisson's equa- 
tion that the potential change across a space-charge layer 
with a given amount of charge is larger when the charge 
is spatially spread over more widely. These features are 
transparent for Fang-Howard wave function, which al- 
lows analytic evaluation as I(z)~z — Q(z)[z — 3/6 + 
cxp(-6z)(3 + 2bz + b 2 z 2 /2)/b] and J(z dcpl )~/(oo)=3/6. 

The net shift in E c across the 2DEG (from z<0 to 
z>0) brought about by electrons themselves is given by 

e 2 

SE C = n e /(z dep i), (13) 

which includes the contribution from the tail of "4>(z) in 
the barrier layer. Although V xc (z) also depends on n e , 
it does not contribute to the net shift since it vanishes 
where ip(z) vanishes (see eqs. (7), (8) and an example in 
Fig. 6). Note that /(zdc P i)=0 for an ideal zero-thickness 
2DEG \tp(z)\ 2 =5(z), therefore eq. (13) is the effect re- 
sulting exclusively from the finite thickness. 

It is worth pointing out that the last term in eq. (11) 
is the only "plastic" term in <j)(z) that can be modified 
externally from the front surface. All the other terms 
are fixed once the sample quality (Na) and the back- 
gate setting (An bg ) are fixed. An attempt to introduce 
potential modulation from the surface, e.g., by the grat- 
ing induces perturbing spatial variation in E c above the 
2DEG, namely, in 8E C . Following the change in 5E c (r) at 
the position r in the x-y plane, n e (r), -0(r, z), and hence 
Ei(r) are also slightly altered to fulfill eqs. (6), (9), (10) 
and (13) with the modified 5E C , thereby resulting in the 
spatial variation in E\ . When the modulation amplitude 
is small enough, as is the case in the present study, the 
conversion ratio from SE C to Ei can be evaluated by the 
derivative dEi/d(SE c ) at the original n e . Further, with 
a plausible assumption that the amplitude of the varia- 
tion in SE r to be introduced from the surface does not 
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depend on properties at and below 2DEG, i.e., on n e , 
ndepi, and Aribg, the amplitude Vq seen by electrons will 
be proportional to dEi/d(5E c ). 

In Fig. 5(b), E\ is plotted against SE C given by eq. 
(13). (To plot all of them within a single frame, E^s 
are negatively offset by the value at n e =2.0xl0 15 m~ 2 
for each Anb g , hence the notation AE\.) The plots are 
for the same sets of n e and An bg as in (a) . Applying ex- 
ternal perturbation to introduce modulation corresponds 
to slightly shifting the point (6E c ,AEi) along the curve 
of constant Anb g , because the perturbation is assumed 
not to affect Anb g . Therefore the slope of the constant 
Anbg curve is the dEi/d(SE c ) to be considered here. It 
can be seen in Fig. 5(b) and clearer in Fig. 5(c), which 
plots 26 dEi/d(5E c ) versus n e , that the slope has a trend 
of becoming smaller with increasing n e if Anb g is kept 
constant, and also with increasing (more positive) Anb g . 
It is to this trend that we ascribe the observed behavior 
of Vq. The decrease of dEi/d(SE c ) by increasing Anb g 
is readily interpretable as an effect of the thickness of 
the 2DEG: the negative shift of 5E C with n e is more 
rapid for thicker 2DEG while the decrease in E\ does not 
change very much (see Fig. 4 or Fig. 5(b)), and therefore 
dE\/d{8E c ) is smaller for thicker 2DEG. The decrease of 
dEi/d(SE c ) with increasing n e for constant Anb g origi- 
nates from more subtle competition between the shift in 
SE C and E\ \ both shift downward with increasing n e , and 
the rate of the shift decreases with increasing n e whose 
rate (the rate of "the rate of shift" ) decelerating slightly 
more rapidly for 6E C , leaving the rate of change for 5E C 
relatively larger than that of E\. 

In Fig. 5, the results of the calculations are plotted 
for all possible sets of n e and Anb g - In the experiment, 
n e was varied either by the back gate alone, namely 
n e —n e0 + Anb g with n e0 fixed, or by illumination with 
the back-gate voltage fixed to zero Anb g =0. The plotted 
points corresponding to the two experimental modes are 
highlighted by large and small solid circles, respectively, 
with n e o set to 2.0xl0 15 m~ 2 . Figure 5(a) shows that 
rms thickness increases with increasing n e when varied 
by the back gate, while the trend is reversed when n e 
is varied by illumination. As demonstrated in Fig. 5(c), 
on the other hand, dE\/ d{5E c ) decreases for both meth- 
ods with the rate of the change much larger when n e is 
driven by the back gate. Thus, dEi/d(SE c ) behaves in 
qualitatively the same manner as Vq under the variation 
of n e by the back gate and by illumination. 

To see how Na affects the problem, we calculate 
dEi/d(5E c ) for two other values of Na through the same 
procedure as was done for Fig. 5(c), and plot in Fig. 7 
the dE\/ d{8E c ) normalized by its value at n e =2.0xl0 15 
m~ 2 . The figure shows that the behavior does not de- 
pend very much on Na when n e is varied by illumination, 
while the value of Na has strong impact on the rate of de- 
crease when the back gate is in use; the decrease is more 
rapid for smaller Na, i.e., for cleaner 2DEG. This is read- 
ily interpretable in terms of its effect on the confinement 
potential. As manifest from eq. (11), the confinement po- 
tential has the same profile as long as n depl =nd cp i — Anb g , 
which should be positive for the confinement to be op- 
erative, remains unchanged. Small Na hence small nd op i 




1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 

n (10 15 m- 2 ) 

Fig. 7. The rate of change dE\/ d(5E c ) (normalized by the value 
at n e =2.0xl0 15 m -2 and An bg =0.0xl0 15 m -2 ) plotted as a 
function of n e for different values of Na ■ Solid and open symbols 
correspond to varying n e by the back gate and by illumination, 
respectively. Stars represent the calculation without V xc (z) in eq. 
(6). 

intensifies the relative importance of the variation of the 
Anbg. Large solid circles in Fig. 5(c) shows that the 
rate of the change of dEi/d(SE c ) versus n e , varied by 
the back gate, increases with increasing (more positive) 
Anbg, namely, with increasing thickness of the 2DEG. 
Since decreasing nd ep i is equivalent to increasing Anb g , 
it also adds to the change rate. The qualitative behavior 
of dEi/d(5E c ) under the variation of n e , Na, and Anbg 
is reproduced in a simple analytic calculation using Fang- 
Howard wave function (see Appendix). 

For quantitative comparison with the experiment, the 
plots in Fig. 7 for Na—1-0 and 2. Ox 10 20 m~ 3 are replot- 
ted in Fig. 2(b). 27 As mentioned before, the rationale of 
comparing Vq with dEi/d(5E c ) is given by the assump- 
tion that the amplitude of the periodic modulation in SE C 
is determined solely by the grating on the surface and in- 
dependent of n e , Na, and Anb g ; from the assumption, 
the proportionality of Vq and dE\j 1 d{5E c ) is expected to 
result. (Note, however, that it is only the proportion- 
ality that can be compared with experiment, since we 
currently do not have a method to directly measure the 
amplitude of the modulation of SE C .) Fig. 2(b) demon- 
strates that the behavior of Vq under back gate and illu- 
mination is explained by that of dE\/ d(5E c ) fairly well. 
The agreement is better for Na—1-0x 10 20 m~ 3 , imply- 
ing that the density of residual impurity in our 2DEG 
wafer may be closer to this value (or smaller). However, 
since we have used simplified model for the confinement 
potential, and have neglected many complications, e.g., 
possible screening by the Si-doped layer, 28 it seems going 
too far to use this comparison for quantitatively accurate 
determination of Na- 

Finally, we discuss the effect of the exchange and cor- 
relation term V xc (z). In Fig. 7 we also plot the result 
of calculation deliberately omitting the term V xc (z) from 
eq. (6). The plot, along with Fig. 2(b), illustrates that ex- 
periments are better described with inclusion of V xc (z). 
The predominant effect of V xc {z) is twofold: it reduces the 
thickness of the 2DEG hence the downward shift of SE C , 
and also reduces E\ by partially alleviating the penalty 
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in energy due to the Hartree term. It turns out by exam- 
ining the process of the calculation that the decrease in 
Ei by V xc (z) dominates the difference with and without 
the term. The effect is more pronounced for larger n e , let- 
ting the decrease of dEx/d(5E c ) with n e more rapid and 
closer to the experiment. This underlines the importance 
of the exchange-correlation effect, albeit the smallness of 
its magnitude, in understanding the behavior of the po- 
tential modulation in the 2DEG. 

5. Conclusions 

We have reported our experimental result that the am- 
plitude Vo of potential modulation decreases with in- 
creasing electron density n e , varied either by the back 
gate or by illumination. Contrary to the intuition, the 
decrease is much more rapid when the back gate is em- 
ployed. We have ascribed the result to the modification in 
the conversion rate dEi/d(5E c ) for the perturbing mod- 
ulation with its source at the surface (or more generally, 
above the 2DEG) to be transmitted to the subband level 
in the confinement potential V(z). The origin of the sub- 
stantial change in dEi/d(6E c ) can be traced back to the 
high sensitivity, particularly notable in a SH 2DEG, of 
the profile of V(z) hence of the envelope function ip(z) to 
the change in the electron density or the back-gate set- 
ting. Therefore we expect much smaller effect in a 2DEG 
formed in a single quantum well. The present result pro- 
vides a prescription for making the modulation ampli- 
tude large in SH 2DEG; that is, it is advantageous to 
make the 2DEG as thin as possible. This is especially 
important for a small period ULSL for which the ampli- 
tude becomes inevitably small. 

Although we have in the present paper confined our 
interest to ULSL, for which the measurement of the mod- 
ulation amplitude Vo is possible, the decrease of Vq with 
increasing n e will also take place likewise in other low- 
dimensional electron systems based on SH 2DEG with 
potential modulation. Moreover, even in plain SH 2DEGs 
without any artificial modulation, the random potential 
landscape seen by electrons, since its main source is ion- 
ized donors located above the 2DEG, will diminish its 
amplitude with increasing n e . The effect should be born 
in mind in interpreting the experiment with variation of 
n e , especially when n e is varied by the back gate. 15,29,30 
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Appendix: Fang-Howard approximation 

The Fang-Howard approximation, although not 
quantitatively quite accurate, provides analytic 
formula, which will be useful for the qualitative 
understanding of the phenomenon. Here, we calcu- 
late dEi/d(SE c ) using Fang-Howard wave function 
■0fh(z). The lowest subband energy is given by the 
sum of the expectation values of the kinetic en- 
ergy and the potential energy, Ei=(T)+(V), with 
(T) = (^ FH (z)|(-ft 2 /2m*)(d 2 /^ 2 )I^FH^))=^& 2 /8m*. 



In evaluating potential energy, we neglect V xc (z) for 
simplicity. Noting that Vfh(z) vanishes at z<0, we 
obtain, using cq. (11), 

(V) = (V>FH(2)|(-e)<K2)|# H (2)) 



thus, 



Ei = Ry* 



n dep i - An bg - 5n e /16 



* 2 J,2 

4 a B b 



+24wa* B 



^de P i-5n e /16 2N A \ 
b b 2 J 



with n depl =nd C pi — Anb g . As mentioned before, 
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Making use of the relation 16,31 &=[487r(n dopl + 
lln e /32)/a| 3 ] 1 / 3 , and being reminded that N A hence 
ndepi, and Anb g are fixed, b can be made the only free 
parameter by eliminating n e , and therefore, 



dE-L 



dE 1 /db 



d(8E c ) d(8E c )/db 



(3a* B /16Tr)b 3 + 21n' - UN A /b 



5 

16 



(4a B /3n)b* + 32n dcpl 

(33/320)n e - (22/15)A A /6 



"depl 



Mepl 



(ll/48)n e 



(A-4) 

In the last equality, the parameter b is replaced back 
again to n e , except for the term including Na, which 
is negligibly (roughly three orders of magnitude) small 
compared with other terms. It is easy to verify from eq. 
(A-4) that dEi/d(8E c ) decreases with increasing n e when 
(i) n e o=n e — Anbg is fixed (simulating the back-gate con- 
trol) and (ii) n dopl is kept constant (simulating the illu- 
mination), and the rate of change is larger for the former. 
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